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Abstract 

We consider the orthogonality graph fi(ra) with 2 n vertices correspond- 
ing to the vectors {0, 1}", two vertices adjacent if and only if the Hamming 
distance between them is n/2. We show that, for n = 16, the stability 
number of fi(n) is a(f2(16)) = 2304, thus proving a conjecture by Gal- 
liard 0. The main tool we employ is a recent semidefinite programming 
relaxation for minimal distance binary codes due to Schrijver |17f . 

As well, we give a general condition for Delsarte bound on the (co)cli- 
ques in graphs of relations of association schemes to coincide with the 
ratio bound, and use it to show that for fl(n) the latter two bounds are 
equal to 2 n /n. 

Keywords: Semidefinite programming, minimal distance codes, stability num- 
ber, orthogonality graph, Hamming association scheme, Delsarte bound. 
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1 Introduction 

The graph Q(n) and its properties 

Let 0(n) be the graph on 2" vertices corresponding to the vectors {0, 1}™, such 
that two vertices are adjacent if and only if the Hamming distance between 
them is n/2. Note that Q(n) is fc-regular, where k — (P ). 

It is known that f2(n) is bipartite if n — 2 mod 4, and empty if n is odd. We 
will therefore assume throughout that n is a multiple of 4. The graph owns its 
name to another description, in terms of ±l-vectors. Then the orthogonality of 
vectors corresponds to the Hamming distance n/2. 
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Moreover, f2(n) consists of two isomorphic connected components, Qo(n), re- 
spectively f2i(n), containing all the vertices of even, respectively odd, Hamming 
weight. For a detailed discussion of the properties of f2(n), see Godsil [5], the 
PhD thesis of Newman i an d EH • 
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Figure 1: The connected component of f2(4) corresponding to vertices of even 
Hamming weight. 

In this note we study upper bounds on the stability number a(0(n)). 
Galliard |7j pointed out the following way of constructing maximal stable sets 
in fl(n). Consider the component O e (n) of fi(n), for 1 — e = j mod 2, and take 
all vertices of Hamming weight e, e + 2, . . . , e + 21, . . . , n/A — 1. Obviously, these 
vertices form a stable set S of Q(n) of size 




We can double the size of S by adding the bit-wise complements of the vertices 
in 5, and double it again by taking the union with the corresponding stable set 
in f2i_ e (n). Thus we find that 

L"/8J , s 

a(fi(n))>4£ L_J :=a(n) (1) 

i— e ^ ' 

For n = 16 this evaluates to a(fl(n)) > 2304. Galliard et al [5] could show that 
a(0(16)) < 3912. In this note we will show that, in fact, a(f2(16)) = 2304. This 
was conjectured by Galliard 0, and Newman ^1] has recently conjectured that 
the value Q actually equals a(fl(n)) whenever n is a multiple of 4. 
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A quantum information game 

One motivation for studying the graph Q(n) comes from quantum information 
theory. Consider the following game from |S]. 

Let r > 1 and n = 2 r . Two players, A and B, are asked the questions xa and 
xb, coded as n-bit strings satisfying 

d H {xA,x B ) e |o, inj 

where dn denotes the Hamming distance. A and B win the game if they give 
answers da and ys, coded as binary strings of length r such that 



VA = Vb 



xa = xb- 



A and B are not allowed any communication (except a priori deliberation). 
It is known that A and B can always win the game if their r output bits are 
maximally entangled quantum bits [2] (see also |14|1. 

For classical bits, it was shown by Galliard et al [5] that the game cannot always 
be won if r = 4. The authors proved this by pointing out that whether or not 
the game can always be won is equivalent to the question 

X(fi(n)) < n? 

Indeed, if \(£l(n)) < n then A and B may color f2(n) a priori using n colors. 
The questions xa and xb may then be viewed as two vertices of f2(n), and the A 
and B may answer their respective questions by giving the color of the vertices 
xa and xb respectively, coded as binary strings of length log 2 n = r. 
Galliard et al. [H] showed that x(^(16)) > 16, i.e. that the game cannot be won 
for n = 16. They proved this by showing that a(f2(16)) < 3912 which implies 



X(Q(16)) > 



>16 



a(n(16)) 



> 



>16 



3912 



17. 



In this note we sharpen their bound by showing that a(fi(16)) = 2304, which 
implies x(^(16)) > 29. 

Our main tool will be a semidcfinitc programming bound on a(Q(n)) that is 
due to Schrijver |17j . where it is formulated for minimal distance binary codes. 



2 Upper bounds on a(Q(n)) 

In this section we give a review of known upper bounds on a(fl(n)) and their 
relationship. 

2.1 The ratio bound 

The following discussion is condensed from Godsil jH]- 
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Theorem 1. Let G = (V, E) be a fc-regular graph with adjacency matrix A(G), 
and let \ m i n (A(G)) denote the smallest eigenvalue of A(G). Then 

I VI 

a(G) < 1 1 . (2) 

1 - A min (A(G)) 

This bound is called the ratio bound, and was first derived by Dclsarte @] for 

graphs in association schemes (see Sect. 12.21 for more on the latter). 

Recall that Q(n) is /c-regular with k = (I 1 )- Ignoring multiplicities, the spec- 

\ 2 n/ 

trum of fi(n) is given by 

25" 

A m = , ,, (to - l)(m - 3) ■ ■ ■ (to - n + 1) (m = l,...,n). (3) 
(2 n ) ! 

The minimum is reached at to = 2, and we get 

25" ( " ) 

A min (^(f2(n)) = _(i)(-i)(- 3 ) • • • (-n + 3) = (4) 
(i?i)! n — 1 

The ratio bound therefore becomes 

on 

a (0(n)) < _. (5) 
n 

This is the best known upper bound on a(f2(n)), but it is known that this bound 
is not tight: Frankl and Rodl [Bj showed that there exists some e > such that 
a(£l(n)) < (2 — e)". For specific (small) values of n one can improve on the 
bound 10), as we will show for n < 32. 

2.2 The Delsarte bound and 

Here we are going to use more linear algebra that naturally arise around fl(n). 
We recall the following definitions, cf. e.g. Bannai and Ito P^. 

Association schemes. An association scheme A is a commutative subalgebra 
of the full v x w-matrix algebra with a distinguished basis (Aq = I, A±, . . . , A n ) of 
0-1 matrices, with an extra property that J^i Ai equals the all-ones matrix. One 
often views Aj, j > 1, as the adjacency matrix of a graph on v vertices; Aj is 
often referred to as the j-th relation of A. As the A/'s commute, they have n+ 1 
common cigcnspaccs Vi. Then A is isomorphic, as an algebra, to the algebra 
of diagonal matrices diag(Po.j> ■ ■ ■ ,-Pnj), where denotes the eigenvalue of 
Aj on Vi. The matrix P = (Pij) is called first eigenvalue matrix of A. The 
set of A,-'s is closed under taking transpositions: for each < j < n there 
exists j' so that Aj = Aj,. In particular, = Pij' ■ An association scheme 
with all Aj symmetric is called symmetric, and here we shall consider such 
schemes only. There is a matrix Q (called second eigenvalue matrix) satisfying 
PQ = QP = vl. In what follows it is assumed (as is customary in the literature) 
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that the eigenspace Vq corresponds to the eigenvector (1, . . . , 1); then the 0-th 
row of P consists of the degrees Vj of the graphs Aj . It is remarkable that the 
0-th row of Q consists of dimensions of Vi . 

Let denote the Schrijver i?'-function |16j : 

tf'(G) = max{Tr(JX) : Tr (AX) = 0, Tr(X) = 1, X h 0, X > 0} . 

For any graph G one has a(G) < i?'(G). Moreover, $'(G) is smaller than or 
equal to the ratio bound for regular graphs, as noted by Godsil Sect. 3.7]. 
For graphs with adjacency matrices of the form J2j^M Aj> wrt h Aic{l,..,,n} 
and A,'s from the 0-1 basis of an association scheme A, the bound coincides, 
as was proved by Schrijver JS| , with the following bound due to Delsarte j3] 0] 

max l T u> subject to w > 0, Q T u> > 0, wo = 1, Wj = for j e M, (6) 

where Q is the second eigenvalue matrix of A. 

The bound © is often stated for (and was originally developed for) bounding 
the maximal size of a g-ary code of length n and minimal distance d; then 
the association scheme A becomes the Hamming distance association scheme 
H(n,q) and M. = {1, . . . , d — 1}. The relations of H(n,q) can be viewed as 
graphs on the vertex set of n-strings on {0, . . . , q— 1}: the j-th graph of H (n, q) 
is given by 

1 if d„{X,Y)=j 
otherwise 



(AO 



For -ff (n, g) the first and the second eigenvalue matrices P and Q coincide, and 
are given by = Kj(i), where Kk is the Krawtchouk polynomial 



k 

71 — X 



K k (x) :=]T(-l)^-l) fc 

3=0 



k-j 



For £l(n), the bound (JSJ) is as above with A = H(n, 2) and Ai = {§}• Newman 
[H] has shown computationally that $'(fi(n)) = 2 n /n if n < 64, i.e. the ratio 
and bounds coincide for fi(n) if n < 64. We show that it is the case for all 
n, as an easy consequence of the following. 

Proposition 1. Let A be an association scheme with the 0-1 basis (Aq, . . . , A n ) 
and eigenvalue matrices P and Q. Let A r have the least eigenvalue r = Pi r and 
assume 

V r Pii > ViT, < i < n. 

Then the Delsarte bound JSJl, with Ai = {?•}, and the ratio bound for A r 
coincide. 

Proof. Let Pj denote the j-th row of P. 

As we already mentioned, the bound J5J for regular graphs always majorates 
(©■ Thus it suffices to present a feasible vector for the LP in (JBJl that gives the 
objective value the same as |j5J. 
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We claim that 



pT , u r p T 



is such a vector. It is straightforward to check that ciq = 1 and a r = 0, as 
required. By the assumption of the proposition, a > 0. As PQ = vl, any 
nonnegative linear combination z of the rows of P satisfies Q T z T > 0. As a T is 
such a combination, we obtain Q T a > 0. 

Finally, to compute l T a, note that 1 t Pq = v and \ T Pj = 0. □ 
Corollary 1. The bounds (JSJ) and J5J) coincide for f2(n). 

Proof. We apply Propositionnto A = H(n, 2) and r = Then the eigenvalues 
of A r = Q(n) given in © comprise the r-th column on P, in particular the least 
eigenvalue r equals P2. r , by J3J above. The assumption of the proposition 
translates into 1 



(n \.. 2f+ 2 ( ? i -2)!(n -!)!!(§ -i) 



2 



^ ( 2) - ^ j A- f (2) = ,, ( » ) ; (n _; ) ," ^ 

as claimed. □ 



2.3 Schrijver's improved SDP-based bound 

Recently, Schrijver J7| has suggested a new SDP-based bound for minimal dis- 
tance codes, that is at least as good as the bound, and still of size polynomial 
in n. It is given as the optimal value of a scmidcfinitc programming (SDP) 
problem. 

In order to introduce this bound (as applied to a(Q(n))) we require some nota- 
tion. 

For i, j, t € {0, 1, ... , n}, and X, Y S {0, 1}" define the matrices 



(Mf 



1 if\X\=i,\Y\=j,d H (X,Y)=i + j-2t 
l >o>x.Y 1 otherwise 



1 f.O 



The upper bound is given as the optimal value of the following semidefinite 
program: 

a(n) := max^ 
i=o 

subject to 

r° - 1 
x o,o — 1 

< x\ 3 < x° for all i, j, t G {0, . . . , n} 

x\ j = x\, j, if {i 1 , j' , i' + j' — 2t'} is a permutation of {i, j, i + j — 2t} 
j 



x\„ = Oif {i,j,i + j-2t}na n }^<D, 



1 Here m\\ = m(m — 2)(m — 4) ... , the double factorial. 



G 



as well as 

The matrices M* • are of order 2™ and therefore too large to compute with in 
general. Schrijver pointed out that these matrices form a basis of the Tcrwilligcr 
algebra of the Hamming scheme, and worked out the details for computing the 
irreducible block diagonalization of this (non-commutative) matrix algebra of 
dimension 0(n 3 ). 

Thus, analogously to the -d'-case, the constraint J2i j t x \ j^i j — ® * s re pl ace d 

by 

i,j,t 

where Q is an orthogonal matrix that gives the irreducible block diagonalization. 
For details the reader is referred to Schrijver ^7j. Since SDP solvers can exploit 
block diagonal structure, this reduces the sizes of the matrices in question to 
the extent that computation is possible in the range n < 32. 



2.4 Laurent's improvement 

In Laurent ^3] one finds a study placing the relaxation ^7] into the framework 
of moment sequences of I12j . This study also explains the relationship with 
known lift-and-project methods for obtaining hierarchies of upper bounds on 
a(G). 

Moreover, Laurent |13| suggests a refinement of the Schrijver relaxation that 
takes the following form: 

l+{n) := max2™XQ 

subject to 

< x\j < x° fi for all i,j,te {0, . . . , n} 
x\ j = x\, j, if {%', j' ', i' + j' — 2t'} is a permutation of {i,j, i + j — 2t} 
xt tj = Oif{i,j,i + i-2t}n{|n}^0, 



as well as 



and 



l-x° n 



0.0 



c 



T 



C ^i,j,t \ X i+j~2t,0 X i,j) 

where c := J2™ =0 { x o,o ~ x o,i) Xi, and \% is defined by 



h0, 



{Xi)x ■-- 



1 \i\X\=i 
else. 



This SDP problem may be block-diagonalised as before to obtain an SDP of 
size 0(n 3 ). 
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3 Computational results 



To summarize, the bounds we have mentioned satisfy: 

a(n) < a{n{n)) < l + (n) < a(n) < tf'(fi(n)) = 2 n /n. 

In Tabled we show the numerical values for a(n) and l+(n) that were obtained 
using the SDP solver SeDuMi by Sturm with Matlab 7 on a Pentium IV 
machine with 1GB of memory. Maltab routines that we have written to generate 
the corresponding SeDuMi input are available online 15 . 



n 


a(n) 


l + (n) 


a(n) 


#'(fi(7i)) = [2 n /n\ 


16 


2304 


2304 


2304 


4096 


20 


20, 144 


20,166.62 


20,166.98 


52,428 


24 


178,208 


183,373 


184, 194 


699,050 


28 


1,590,376 


1,848,580 


1,883,009 


9,586,980 


32 


14, 288, 896 


21,103,609 


21,723,404 


134,217, 728 



Table 1: Lower and upper bounds on a(Q(n)). 

Note that the lower and upper bounds coincide for n = 16, proving that 
a(f2(16)) — 2304. The best previously known upper bound, obtained by an 
ad hoc method, was a(fi(16)) < 3912 ||. 
The value a(20) = 20, 166.98 implies that 

a(fi(20)) £ {20144,20148,20152,20156,20160,20164} 

since a(fl(n)) is always a multiple of 4. Another implication is that n = 20 is 
the smallest value of n where the upper bounds a(n) and l+(n) are not tight. 
It is worth noticing that the Schrijver and Laurent bounds (ot(n) and l+(n) 
respectively) give relatively big improvements over the Delsarte bound — . This 
is in contrast to the relatively small improvements that these bounds give for 
binary codes, cf. [171 113| . We also note that these relaxations are numerically 
ill-conditioned for n > 24. This makes it difficult to solve the corresponding 
SDP problems to high accuracy. The recent study by De Klerk, Pasechnik, 
and Schrijver [2] suggests a different way to solve such SDP problems, leading 
to larger SDP instances, but which may avoid the numerical ill-conditioning 
caused by performing the irreducible block factorization. 
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